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1.  Introduction  : 

II  y a encore  queiques  annees,  le  developpement  des  methodes  et  des  outils  de  calcul  en 
aeroelasticite  constituait  chez  Dassault  aviation  un  challenge  important.  L’outil  de  calcul 
CATIA-ELFINI-AEROELASTICITY  connaissait  a cette  epoque  un  essort  important,  soutenu 
par  des  budgets  et  des  equipes  de  developpement  consequents,  ainsi  que  par  des  besoins 
grandissants  dans  le  domaine  du  design  et  de  la  justification  d’avions  modernes  tels  que 
RAFALE  ou  FALCON  2000. 

Aujourd’hui,  les  priorites  changent  et  visent  principalement  a parvenir  a maitriser 
completement  les  outils  complexes  disponibles,  plutot  que  de  rechercher  de  nouvelles 
methodes.  Des  strategies  industrielles  de  calculs  et  d’essais  sont  mises  au  point  autours  des 
outils  existants.  Elies  ont  pour  but  de  permettre  a des  ingenieurs  de  calculs  specialises,  rnais 
aussi  a des  ingenieurs  de  conception  dans  les  bureaux  de  dessin  ou  a des  ingenieurs  d’essais 
en  vol,  de  disposer  de  modeles  fiables  au  plus  tot  dans  le  cours  du  projet  et  a reduire  le 
nombre  de  configurations  a tester,  ainsi  que  le  nombre  d’essais  en  vol.  L’objectif  clairement 
affiche  est  la  maitrise  des  risques  et  des  couts  en  aeroelasticite  dans  le  cadre  du  projet  avion. 
Dans  les  sections  qui  suivent,  nous  illustrons  les  points  les  plus  importants  de  ces  travaux 
chez  Dassault  Aviation.  Nous  mettons  ainsi  en  relief : 

• L’ ensemble  des  methodes  de  calcul  aujourd’hui  disponibles  en  aeroelasticite. 

• La  mise  ne  oeuvre  de  ces  techniques  dans  le  cadre  de  la  maitrise  des  risques  et  de  la 
reduction  des  couts. 

• A plus  long  terme,  les  techniques  d’ optimisation  indispensable  pour  permettre  a 
l’avionneur  de  prendre  en  compte  rationnellement  dans  ces  plans  de  qualification  et  dans 
ses  dossiers  de  justification  les  problemes  lies  a la  multiplicity  et  a la  variability  des 
configurations. 

• L’ utilisation  systematique  de  methodes  de  calibration  des  procedures  de  calculs  et  des 
modeles  numeriques  sur  des  resultats  d’essais  reels. 

Plus  concretement,  nous  donnons  des  elements  sur  la  qualite  que  l’on  peut  attendre 
aujourd’hui  des  previsions,  ainsi  que  les  gains  esperes  en  matiere  de  reduction  du  nombre 
d’essais  en  vol. 


2.  Outils  et  methodes  disponibles  en  Aeroelasticite 

2. 1 Etat  de  I’art  en  aeroealsticite  lineaire 

Sur  le  plan  theorique,  l’etat  de  Part  en  aeroelasticite  lineaire  est  encore  globalement  fonde  sur 
la  methodologie  lineaire  classique,  dont  les  principes  etaient  deja  bien  etablis  dans  les  annees 
60.  Le  detail  des  methodes,  les  strategies  de  calcul  et  les  techniques  de  calibration  associees 
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ont  ete  largement  decrits  dans  la  Ref  [1],  Nous  en  resumons  les  principales  caracteristiques 
dans  cette  section. 

Les  modeles  mathematiques  intervenant  dans  la  chaine  de  calcul  aeroelastique  sont  decrits  ci- 
dessous  : 
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Figure  LAeroelastic  analys  process 

On  utilise  un  modele  elastique  elements  finis  (FEM)  de  F avion  complet  comprenant  de 
l’ordre  de  10000  a 100000  degres  de  liberte.  Pour  la  modelisation  aerodynamique  une 
modelisation  en  potentiel  lineaire  par  la  methode  des  singularity  est  etablie  dans  le  domaine 
frequence.  On  pallie  les  imprecisions  de  ce  modele  en  recalant  certains  termes  par  des 
resultats  de  calculs  CFD  ou  en  les  recalant  sur  des  resultats  d’essais  en  soufflerie  ou  d’essais 
en  vol. 

Le  principe  est  alors  d’effectuer  les  operations  lourdes  sur  le  modele  elements  finis  et  sur  le 
modele  aerodynamique  independamment  de  toutes  les  configurations  de  vol  et  de  toutes  les 
repartitions  de  masses  a considered  Ces  operations  lourdes  sont  les  suivantes  : 

• Le  calcul  des  deplacements  du  modele  elements  finis  pour  quelques  centaines  de 
chargements  de  base.  Ces  chargements  de  base  correspondent  a des  charges  unitaires  de 
pression  ou  d’inertie  par  zones.  Tout  chargement  statique  ou  dynamique  peut  etre 
represente  par  combinaison  lineaire  des  chargements  de  base. 

• La  construction  du  modele  aerodynamique  exprimant  les  charges  aerodynamiques 
stationnaires  et  instationnaires  dans  la  base  de  charges  precedentes,  en  fonction  du 
mouvement  de  1’ avion.  Ce  mouvement  est  defini  par  une  combinaison  de  formes  de  base 
aerodynamiques  independantes  des  degres  de  liberte  du  modele  elastique. 

• L’ extraction,  a partir  de  ces  modeles  elastiques  et  aerodynamiques  encombrants  des 
‘operateurs  de  base  condenses’  nontenant  les  seules  donnees  intervenant  dans  le  couplage 
aeroelastique  : matrices  de  rigidite  et  de  masses  reduites  dans  la  base  de  charge,  operateurs 
de  lissage  des  deformees  elements  finis  de  base  par  les  forces  aerodynamiques  de  base, 
‘indicateurs  de  severite’  des  charges  de  base. 

A partir  de  ces  operateurs  de  base,  les  calculs  qui  suivent  sont  tres  rapides  : 

• L’elaboration  des  coefficients  des  equations  de  l’aeroelasticite  pour  les  configurations  de 
vol  et  de  masse  considerees. 

• Les  analyses  de  stabilite  (divergence  et  flutter). 

• Le  calcul  du  mouvement  souple  de  F avion  en  quasi-statique,  ou  en  dynamique  (incluant  la 
prevision  des  fonctions  de  transfert  pour  les  etudes  de  stabilite  CDVE)  avec  le  suivi  des 
indicateurs  de  severite. 

• La  selection  des  cas  de  charges  dimensionnants. 
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Figure  2:  Linear  aeroelasticity  typical  outputs 


2.2  Prise  en  compte  des  non-linearites 

Afin  d’illustrer  les  methodes  aujourd’hui  disponibles  et  industrielles  dans  ce  domaine,  nous  ne 
donnons  ici  que  quelques  exemples  typiques  de  traitement  des  non-linearites  aeroelastique. 


• Non-linearites  aerodynamiques  : 

En  pratique,  que  ce  soit  en  stationnaire  ou  en  instationnaire,  la  methode  lourde  consiste 
aujourd’hui  a iterer  avec  une  analyse  d’aeroelasticite  lineaire  recalee  a chaque  iteration  sur 
des  calculs  CFD  (full-potentiel,  Euler,  ...)•  Les  efforts  aerodynamiques  provoques  par  les 
deformations  souples  sont  alors  calcules  par  la  methode  classique  des  singularity  lineaires  et 
servent  de  pre-conditionneur.  La  convergence  est  alors  acquise  en  quelques  iterations  comme 
le  montre  la  figure  3 (calcul  du  Cz  d’un  avion  FALCON  equilibre  a lg  - M=0.85  - Couplage 
CATIA/ELFINI/AEROELASTICITY  - EULER/EUGENIE)  : 


CSM  : Catia-Elfini- Aeroelasticity  Euler  -CFD : Eugenie 


Pitch  moment  (m3)  for  a unitary  angle  of  attack  maneuver 
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Figure  3:  Example  of  iterative  CSM-CFD  coupling  process 


Meme  si  1’  augmentation  des  performances  des  moyens  de  simulation  numeriques  permet 
d’envisager  dans  Favenir  F utilisation  plus  courante  d’outils  CFD  en  aeroelasticite,  le  nombre 
de  cas  de  calculs  a balayer  pour  couvrir  l'ensemble  du  domaine  de  vol  et  1'ensemble  des 
configurations  structurales  restent  trop  importants  pour  pouvoir  envisager  l'utilisation  massive 
de  ces  outils  a moyen  terme.  On  utilise  done  plutot,  lorsque  cela  est  possible,  les  outils  de 
CFD  linearises  qui  permettent  d’integrer  les  resultats  de  calcul  CFD  dans  la  structure  de  la 
base  de  donnees  aeroelastique  lineaire  classique  (Cf.  § 2.1),  ce  qui  permet  d'en  conserver  les 
performances  (cout,  stabilite,  precision).  A titre  d’illustration,  la  planche  suivante  donne 
l’exemple  de  forces  generalisees  instationnaires  calculees  sur  les  premiers  modes  de 
flexion/torsion  voilure  d’un  FALCON  2000,  par  couplage  de  l’outil  CATIA-ELFINI  avec 


60-4 


l’outil  EUGELIN  (outil  CFD  Euler  linearise  - Cf.  Ref  [2]).  Elies  sont  comparees  aux  memes 
forces  generalisees  calculees  par  la  methode  lineaire  classique  des  singularity. 


Generalized  unsteady  aerodynamic  forces  Classical  flutter  analysis  (PK  method) 

CSM  (ELF INI)  / Linearized  CFD  (EUGENIE) 

Figure  4:  Example  of  generalized  unsteady  aerodynamic  forces  computed  using  linearized  CDD 


• Non-linearites  mecaniques  : 

Dans  ce  domaine,  les  types  de  non-linearites  les  plus  souvent  rencontrees  sur  avion  sont 
principalement  dues  : 

> Aux  effets  des  grandes  rotations  de  1’ avion  (mecanique  du  vol), 

> aux  non-linearites  de  contact  liees  a la  presence  de  jeux  mecaniques  (rails, 
systemes  Fail-Safe,  timonerie  de  servo-commandes,  ...), 

> aux  couplages  de  F avion  avec  des  systemes  mecaniques  a comportement 
non-lineaires  (train  d’ atterrissage,  servo-commandes,  ejecteurs  d’emports, 
...). 

Sur  avion  militaire,  l’exemple  typique  illustrant  la  prise  en  compte  de  ces  phenomenes  est  la 
modelisation  des  phases  d’ ejection  en  vol  de  bidons  sous  voilure  de  RAFALE  (Cf.  figure  5) 
ou  la  modelisation  des  phases  de  tirs  de  missiles  sur  rail  au  point  extreme  de  voilure  du 
RAFALE,  pendant  des  manoeuvres  complexes  de  F avion  (Cf.  figure  6). 


Figure  5:  In-flight  tank  ejection  simulation 
using  CATIA-ELFINI-AEROELASTICITY 


Figure  6:  In-flight  wing-tip  missile  launch  simulation 
using  CATIA-ELFINI-AEROELASTICITY 


Dans  ces  deux  cas,  F outil  de  calcul  doit  pouvoir  rendre  compte  avec  une  bonne  precision  du 
couplage  simultane  des  3 types  de  non-linearites  precedentes  dans  la  simulation. 
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II  est  a noter  que  ces  analyses  non-lineaires  necessitent  des  integrations  dans  le  domaine 
temporel  avec  formulation  prealable  des  forces  aerodynamiques  instationnaires  dans  ce 
domaine.  Les  techniques  utilisees  pour  cela  sont  issues  des  travaux  de  Karpel  (Cf.  Ref.  [3]) 
sur  la  rationalisation  des  forces  aerodynamiques  instationnaires  (Cf.  figure  7). 


Unsteady  generalized  forces 
(real  and  imag  part) 


Frequency  Doublet-lattice  (Using  Karpel  rationalized  afro)  (Using  doublet  lattice  atro) 

Figure  7 : Rationalization  of  unsteady  aerodynamic  forces  using  Karpel  method 


2.3  Recherche  des  configurations  les  plus  critiques 

Les  moyens  d’ analyse  de  l’aeroelasticite  que  nous  venons  de  presenter  dans  les  sections 
precedentes  sont  extremement  utiles  pour  le  calcul  d’une  configuration  particuliere  de  1’ avion. 
Cependant,  l’activite  de  calcul  aeroelastique  se  traduit  par  un  nombre  important  (voire  infini) 
de  configurations  structurales  et  massiques  a etudier. 

En  general,  le  balayage  des  cas  de  calcul  s’effectue  sur  un  choix  a priori  des  configurations  les 
plus  critiques.  La  pertinence  de  ce  choix  repose  sur  l’experience  du  calculateur  et  sur  sa 
comprehension  des  phenomenes  physiques. 

Dans  le  cas  des  calculs  de  flutter,  nous  avons  ete  amene  a developper  un  outil  d’ optimisation 
pour  la  recherche  des  configurations  les  plus  critiques  (configurations  dites  ‘les  pires’),  afin 
de  guider  le  raisonnement  du  calculateur,  limiter  les  risques  d’oublis  et  appuyer  la 
demonstration  qui  conduit  au  choix  de  ces  configurations.  Cet  outil  est  egalement  utilise  pour 
permettre  de  tester  la  sensibilite  des  analyses  de  flutter  vis  a vis  des  parametres  incertains  des 
modelisations  ( parametres  ‘flous’  des  modelisations : position  du  reservoir  dans  les 
reservoirs,  rigidite  des  interfaces,  interactions  aerodynamiques  complexes,  ...). 

Les  principes  theoriques  de  cet  outil  ont  ete  abondamment  decrits  dans  la  reference  [4], 


3.  Exemple  de  retombees  des  etudes  aeroelastiques  sur  la  maitrise 
des  risques  et  des  couts  : etudes  dites  ‘de  transposition’ 

Les  etudes  dites  ‘de  transposition’,  mises  en  place  chez  Dassault  Aviation  ces  dernieres 
annees  sur  les  programmes  MIRAGE  2000  et  RAFALE,  concernent  la  stabilite  aeroelastique 
(flutter)  et  aeroservoelastique  (accrochage  du  systeme  de  commande  de  vol  sur  la  structure) 
de  1’ avion.  Elies  ont  pour  objectifs  : 

• De  predire  par  calcul,  avec  un  niveau  de  precision  suffisant,  le  flutter  et  les  fonctions  de 
transfert  souples  CDVE1,  pour  les  configurations  de  base  de  1’ avion. 

• De  verifier  et  de  recaler  eventuellement  ces  previsions  sur  la  base  d’essais  de  vibrations  au 
sol  et/ou  en  vol. 


1 On  appelle  fonction  de  transfert  CDVE  les  fonctions  de  transfert  des  deformations  souples  au  niveau  des 
capteurs  gyrometre  et  accelerometre  du  systeme  de  commande  de  vol  sur  les  braquages  de  gouvernes. 
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• De  transposer  ees  resultats  sans  essais  supplementaires,  pour  toutes  les  versions  d’avions 
(proto  <->  serie,  monoplace  <->  biplace,  Air  version  <->  Marine  version), dans  tout  le 
domaine  de  vol  et  pour  toutes  les  sous-configurations. 

Ceci  suppose  de  disposer  de  modeles  capables  de  predire  les  instability  aeroelastiques  et  les 
fonctions  de  transfert  souples  CDVE,  dans  tout  le  domaine  de  vol,  quelque  soit  la  complexity 
de  la  configuration  et  avec  un  niveau  de  confiance  maximal.  Ceci  suppose  egalement 
d’ avoir  maitrise,  integre  et  valide  des  methodes  de  recalage  experimentales  de  modele 
efficaces  (Cf.  §4). 

Les  retombees  importantes  attendues  de  ces  etudes  concernent : 

• La  verification  avant  vol  des  marges  de  stabilite  de  1’ avion  : retombees  sur  la  maitrise 
des  risques 

• La  definition  des  ouvertures  de  domaine  de  vol  au  ‘juste  assez’  (definition  des 
progressions  des  vols  d’ouverture  suivant  les  marges  aux  calculs)  : retombees  sur  la 

maitrise  des  couts. 

Le  chemin  critique  de  cet  organigramme  est  constitue  par  les  delais  (souvent  trop  courts)  entre 
la  fin  des  essais  au  sol  (dont  le  TO  est  soumis  a l’approvisionnement  des  maquettes  d’emports) 
et  le  debut  des  essais  en  vol  d’ouverture  de  domaine.  Operationellement,  les  plannings 
d’ essais  en  vol  doivent  done  etre  organises  pour  laisser  suffisamment  de  temps  a 
l’aeroelasticien  pour  calibrer  les  modeles  elastiques  sur  la  base  des  essais  de  vibrations  au  sol 
de  l’avion.  L’ exploitation  des  vols  d’ouverture  de  domaine  donnera  lieu,  apres  coup,  a un 
recalage  de  la  partie  aerodynamique  instationnaire  du  modele  aeroelastique.  Ce  modele  sera 
alors  utilise  pour  transposer  les  resultats  d’essais  a tout  le  domaine  de  vol,  aux  sous- 
configurations  et  aux  autres  types  d’avions  (proto  <->  serie,  monoplace  <->  biplace,  Air 
version  <->  Marine  version),  ce  qui  permettra  d’ecrire  les  dossiers  de  justification  de  la 
definition  de  l’avion. 

Afin  d’illustrer  la  precision  que  Ton  peut  attendre  aujourd’hui  des  modeles  dans  ce  domaine, 
la  figure  ci-dessous  donne  les  predictions  de  la  fonction  de  transfert  ‘Pitch/symmetric  control 
surface  deflection’  du  RAFALE  au  sol  en  presence  de  missiles  aux  points  externes  et  medians 
de  voilure  et  de  3 bidons  lourds  sous  voilure  et  fuselage.  Ces  calculs  sont  compares  a des 
resultats  d’essais  de  vibrations  au  sol  dans  2 configurations  : bidons  pleins  et  bidons  vides. 


Full  Tanks  (fuselage  + wing)  :Q/Dm 


Empty  Tanks  (fuselage  + wing) : Q/tlm 


Figure  8:  FCS  transfer  function  prediction  (on  the  ground) 

Coinme  le  montrent  ces  comparaisons,  la  precision  du  calcul  est  tres  bonne.  Le  modele  predit 
bien  en  particulier  le  couplage  des  modes  ‘2  nceuds  fuselage’  et  ‘tangage  missile  au  point 
extreme  de  voilure’,  a 9.5  Hz  et  1 1.5  Hz  lorsque  les  bidons  sont  pleins,  en  un  seul  et  meme  pic 
de  frequence  1 1 hz,  lorsque  les  bidons  se  vident.  Connne  le  montrent  le  modele  et  les  essais, 
ce  couplage  conduit  a une  augmentation  de  la  fonction  de  transfert  de  nature  a diminuer  la 
marge  de  stabilite  des  commandes  de  vol.  Dans  cette  affaire,  la  bonne  prevision  des  modeles 
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nous  a permis,  sur  la  base  des  calculs,  de  reduire  le  nombre  de  vols  d’ouverture  de  domaine 
initialement  prevus  d’environs  50  % (6  vols  realises  au  lieu  des  12  initialement  prevus). 

Autre  exemple  d’ illustration  : la  prediction  de  l’influence  d’un  lest  augmente  en  pointe  avant, 
sur  la  fonction  de  transfert  en  vol  (Mach=0.6)  ‘Pitch/symmetric  control  surface  deflection’  du 
RAFALE  en  configuration  Air/ Air  : 


Figure  9:  FCS  transfer  function  prediction  (on  the  ground) 

Comme  le  montre  ces  courbes,  le  modele  aeroelastique  predit  bien  la  diminution  en  frequence 
du  mode  ‘2  noeuds  fuselage’  due  a 1’ augmentation  du  lest.  Cette  diminution  en  frequence  se 
traduit  a contrario  par  une  augmentation  de  F amplitude  de  ce  pic  en  vol,  qui  peut  etre  source 
d’accrochage  des  commandes  de  vol  sur  la  structure  si  ce  phenomene  n’est  pas  correctement 
maitrise. 


4.  Validations  experimentales  des  procedures  de  calculs  et  des 
modeles 


Comme  il  a ete  vu  dans  les  sections  precedentes,  F utilisation  des  outils  de  predictions 
aeroelastiques  dans  le  cadre  de  la  maitrise  des  risques  et  des  couts  suppose  que  l’on  progresse 
en  parallele  sur  les  techniques  experimentales  pour  la  validation  des  procedures  de  calcul  et 
des  modeles.  Ceci  suppose  egalement  un  niveau  de  maitrise  comparable  des  methodes 
experimentales  au  sol  et  en  vol. 

La  strategic  aujourd’hui  retenue  en  la  matiere  chez  Dassault  Aviation  est  globalement  la 
suivante  : 

• Utilisation  de  maquettes  souples  en  soufflerie  pour  la  mise  au  point  des  outils,  leurs 
validations  et  la  formulation  des  procedures  d’emplois.  Cette  tache  supporte  egalement  le 
choix  des  methodes  et  des  procedures  d’analyse. 

• Utilisation  d’essais  ‘grandeur  nature’  sur  avions  au  sol  et  en  vol  pour  le  recalage  des 
modeles  aeroelastiques  qui  serviront  dans  le  cadre  des  programmes  avions. 

4.2  Maquettes  souples  en  soufflerie 

Le  recours  a des  essais  en  soufflerie  sur  maquettes  souples  nous  parait  capital.  II  permet  de 
mettre  en  correspondance  mesures  et  modelisations,  afin  de  valider  les  procedures  et  les 
methodes  de  calculs  en  aeroelasticite.  II  nous  est  en  effet  apparu  tres  clairement  dans  le  passe 
que  F utilisation  de  resultats  d’essais  reels  permet  de  faire  progresser  et  de  valider  les 
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techniques  de  calcul,  en  devoilant  des  difficultes  non  decelables  dans  le  cas  de  validations 
purement  numeriques. 

Dans  la  mesure  du  possible,  on  s’arrange  pour  que  l’organisation  de  ces  essais  suive  une 
trame  semblable  a celle  de  l’enchainement  conception/calculs/essais  d’un  projet  avion  : 


Design 

Computation 

Manufacturing 

Ground  Static  Test 
Gauge  Calibration 
GVT 


Design 

Computation 

Manufacturing 

Ground  Static  Test 
Gauge  Calibration 
GVT  (Clamped  Tail) 


, N ^Aircraft  development  Study  logic 


FEM  Calibration 
Computation  Update 


FEM  Calibration 
Test  prediction  Update 


Flight  Test  X j j Wind  Tunnel  Test 

Unsteady  Aerodynamic  )j  Unsteady  Aerodynamic 
adjustment  >■  adjustment  (from  WTT) 


Figure  10:  Flexible  Wind  Tunnel  mock-up  study  logic 


Un  exemple  typique  d’ etude  experimentale  menee  dans  ce  cadre  a ete  realisee  en  cooperation 
avec  l’ONERA,  EADS  et  DASSAULT  dans  une  etude  financee  par  le  MOD  Frangais.  II 
s’agissait  de  l’etude  en  soufflerie  (S2  - Modane)  d’une  maquette  souple  d’empennage 
horizontal  d’avion  de  type  AIRBUS  dotee  d’une  gouverne  de  bord  de  fuite,  en  presence  de 
jeux  mecaniques.  Cette  etude  avait  pour  but  de  valider  notre  methodologie  complete  de  calcul 
de  flutter  en  presence  de  jeu  mecanique  intervenant  pour  la  certification  des  timoneries  Fail- 
Safe  des  empennages  de  nos  avions  FALCON.  La  figure  15  donne  un  extrait  des  resultats 
principaux  de  ces  travaux.  Le  detail  complet  de  la  methodologie,  du  montage  et  des  resultats 
d’essais/calculs  peut  etre  trouve  dans  la  reference  [5]. 


Nous  avons  montre  dans  cette  etude  qu’en  procedant  a partir  des  modeles  structuraux  et 
aerodynamiques  recales  classiquement,  la  simulation  restituait  les  resultats  d’essais  a la  fois 
qualitativement  (comportement  amorti,  cycles  limites  et  divergences)  et  quantitativement 
(niveaux  d’ acceleration  des  cycles  limites  provoques  par  les  phenomenes  de  contact  de  l’axe 
de  rotation  de  la  gouverne  sur  les  butees  elastiques). 

4.3  Essais  reels  sol  et  vol  sur  avion 

Ces  essais  sont  utilises  pour  calibrer  et  valider  les  modeles  aeroelastiques  avions  dans  le  cadre 
des  etudes  programmes.  Ils  sont  effectues  pour  une  configuration  de  base  et  quelques 
configurations  identifiees  comme  les  plus  sensibles  : on  admet  alors  qu’apres  calibration  sur 
ces  essais,  le  modele  aeroelastique  puisse  representer  les  autres  configurations. 
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Les  essais  au  sol  sur  avion  (statiques  ou  vibrations  sol)  ont  pour  objectif  de  calibrer  le  modele 
elastique  de  1’ avion.  Les  caracteristiques  modules  identifies  lors  de  ces  essais  peuvent  etre 
utilisees  pour  remplacer  le  modele  dynamique  theorique  dans  les  calculs  de  flutter  et  ainsi 
evaluer  les  differences  de  vitesses  critiques  obtenues  par  rapport  au  calcul  utilisant  le  modele 
elements  finis.  Les  essais  en  vol  (manoeuvres,  vibrations  par  les  gouvemes)  serviront  eux  a 
calibrer  le  modele  aerodynamique  a travers  la  mesure  des  parametres  de  vol,  des  accelerations 
rigides  et  des  reponses  de  la  structure  (jauges  et  accelerations). 


w 

r 


Predictions 


s 


Structural  F.E.  model 
adjustments 


Ground  Tests 


Unsteady  aerodynamic 
model  adjustments 


Flight  Tests 


Figure  12:  Test  versus  model  calibration  strategy 


Les  techniques  associees  de  recalage  des  modeles  elastiques  sur  les  ‘modes  mesures’  et  de  la 
partie  stationnaire  des  charges  aerodynamiques  ont  ete  developpees  et  integrees  dans  notre 
code  CATIA-ELFINI-AEROELASTICITY  dans  les  annees  80.  Elies  sont  aujourd’hui 
utilisees  industriellement  dans  le  cadre  des  programmes  avions  et  ont  fait  1’objet  de  nombreux 
articles  de  notre  part  qui  les  detaillent  (Cf.  Ref.  [6]  et  [7]  par  exemple).  Ces  outils  sont  fondes 
sur  une  technique  mathematique  originale  d’ identification  de  modeles. 

Ces  dernieres  annees,  notre  effort  a essentiellement  porte  sur  les  outils  de  recalage  des  charges 
aerodynamiques  instationnaires,  sur  la  base  des  essais  de  vibrations  en  vol  (Cf.  figure  16). 
Dans  ces  methodes,  nous  utilisons  toujours  la  meme  approche  mathematique  (Cf.  Ref.  [8])  : 

• Les  parametres  d’ajustement  X,  inconnues  du  probleme,  sont  soit  les  forces  generalisees 
instationnaires,  soit  directement  les  composantes  des  champs  de  pression  sur  la  grille 
aerodynamique. 

• La  fonction  de  cout  a minimiser  est  l’ecart  entre  les  valeurs  X recherchees  et  leurs  valeurs 
theorique  L,i,  dans  le  modele. 

• Les  contraintes  a respecter  sont  les  equations  de  restitution,  par  le  modele,  des  frequences 
et  des  amortissements  des  modes  aeroelastiques  mesures  a une  precision  e donnee. 

Afin  d’illustrer  le  fonctionnement  de  cette  methode  dans  le  cadre  du  recalage  des  champs  de 
pressions  sur  la  grille  aerodynamique,  nous  1’appliquons  au  cas  du  calcul  en  flutter  de  la 
maquette  en  soufflerie  souple  decrite  en  Ref.  [5].  Le  diagramme  a gauche  de  la  figure  18 
donne  la  prevision  de  calcul  comparee  aux  mesures  en  soufflerie.  On  constate  que  la  vitesse 
de  flutter  est  surestimee  par  le  calcul  par  rapport  a la  mesure.  Le  diagramme  de  droite  de  la 
figure  18  donne  la  courbe  de  flutter  apres  recalage  des  champs  de  pressions  instationnaires  a 
partir  des  deux  premiers  points  de  frequences  et  amortissements  mesures  a 37800  Pa  et  60000 
Pa  (points  nettement  inferieurs  a la  vitesse  critique).  Comme  le  montre  ce  diagramme,  le 
calcul,  recale  a Laide  de  points  de  mesures  nettement  inferieurs  a la  vitesse  de  flutter,  predit 
desormais  la  bonne  vitesse  critique  a 84500  Pa. 
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Test  / Computation  comparison 
Before  aerodynamic  adjustment 


Frequency  (Hz) 


Test  / Computation  comparison 
After  aerodynamic  adjustment 


Adjustment  2 1 


1.0  2.0 


Figure  13:  Example  of  unsteady  aerodynamic  adjustment  on  WT  tests 


Un  autre  exemple  de  resultat  de  cette  methode  a ete  obtenu  lors  du  recalage  des  forces 
generalisees  instationnaires  du  MIRAGE  FI.  La  configuration  d’etude  est  representee  sur  la 
figure  suivante  : 


Experimental  structural  model  Theoretical  aerodynamic  model 


Figure  14:  Mirage  F 1 : configuration  studied 


A partir  de  cette  configuration,  un  calcul  couplant  les  modes  mesures  au  sol  de  cette 
configuration  a un  modele  theorique  de  forces  aerodynamiques  instationnaires  par  la  methode 
des  doublets  a ete  realise.  La  comparaison  essais/calculs  en  flutter  est  presente  sur  la  figure 
suivante.  Comme  le  montre  la  planche  precedente,  le  calcul  de  flutter  presente  une  vitesse 
critique  nettement  superieure  a celle  des  essais. 
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Figure  15:  Mirage  FI  : comparison  test/prediction 
before  adjustement 


Adjustment  domain 

4 

1 . » © 3 oocuds  voiiure 

2 m X Tangagc  anti.R550 

Point  i 0.032  iu  lieu  dc  0.052 

% 

Figure  16:  Mirage  FI  : comparison  test/prediction 
after  adjustement 


Pour  ameliorer  les  predictions,  ce  calcul  a donne  lieu  a un  recalage  des  forces  generalises 
aerodynamiques  instationnaires  du  modele,  suivant  la  methode  que  nous  venons  de  decrire. 
Pour  cela,  les  mesures  des  frequences  et  des  amortissements  des  trois  premiers  points  de  vol 
(200  m/s,  300  m/s,  400  m/s,  de  vitesses  nettement  inferieures  a la  vitesse  critique)  sont 
utilisees.  Le  calcul  de  flutter  a 1’aide  de  ce  modele  recale  est  presente  ci-dessous.  On  constate 
que  la  vitesse  critique  calculee  a l’aide  du  modele  recale  se  superpose  maintenant  a celle  issue 
des  essais  en  vol. 
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